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Starting from an initial pure quantum state, we present a strategy for reaching a target state 
corresponding to the extremum (maximum or minimum) of a given observable. We show that a 
sequence of pulses of moderate intensity, applied at times when the average of the observable reaches 
its local or global extremum, constitutes a strategy transferable to different control issues. Among 
them, post-pulse molecular alignment and orientation are presented as examples. The robustness 
of such strategies with respect to experimentally relevant parameters is also examined. 
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I. INTRODUCTION 

Recent advances in laser technology have opened up 
large possibilities for the control of quantum processes 
playing a role in a variety of problems encompassing 
chemical reactivity 0, Q, high-order harmonic genera- 
tion Q , or logical gates with applications in quantum in- 
formation 0, M El ■ Several advances have been achieved 
in this context, from the discovery of elementary basic 
mechanisms of strong field induced molecular dynamics, 
to optimal control strategies bringing together sophis- 
ticated numerical algorithms (such as evolutionary al- 
gorithms) 0, Q- Parallel to these works, a systematic 
study of quantum systems in either pure or mixed (sta- 
tistical) states has been undertaken using control theory 
@ . Upper and lower bounds corresponding to kinemati- 
cal constraints have been established for the expectation 
value of arbitrary observables of driven quantum systems 
0- llOj , with a particular attention paid to their dynam- 
ical realizability jTTL . This latter question, which is 
related to their complete controllability, has been solved 
for linear control, i.e., for Hamiltonians depending in a 
linear way on the control functions, of non-dissipative 
finite-level quantum systems 

The aim of this study, which is a continuation of our 
previous paper |13|. is to show how a strategy worked 
out for the laser control of a specific observable can be 
transposed to a generic system. An essential feature of 
the control scheme we propose is to reduce the physical 
Hilbert space to a subspace H^ N \ of finite dimension N. 
A target state |x^)> which maximizes (or minimizes) 
the expectation value of an observable O, is then defined 
in this subspace. The first obvious advantage of the fi- 
nite dimensionality being that |x ) appears, through 
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a variational principle, to be the eigenvector (for the 
non-degenerate case) associated to H^ N ' with the highest 
or lowest eigenvalue of the restriction of the observable 
O. These bounds being kinematic, they only indicate to 
which extent the system can be dynamically controlled 
by application of time-dependent external fields. More- 
over, in the finite-dimensional case, the complete or non- 
complete controllability can be determined and related 
(for a linear control system) to the dimension of the Lie 
algebra generated by the unperturbed Hamiltonian and 
the interaction operator |lll [12( • 

Assuming the dynamical realizability of the target 
state, the remaining question consists in explicitly finding 
a suitable scheme in terms of a unitary evolution opera- 
tor U (t) which, starting from an initial pure state of the 
molecular system, dynamically brings it close to the state 
|x^)- Several control strategies have already been pro- 
posed using, for instance, adiabatic passage techniques 
Il5j . factorizations of unitary operators [l^ or op- 
timal control schemes 0,^3- Here, we pursue another 
approach and we show that a sequence of short pulses of 
moderate intensity can be devised for reaching these tar- 
get states. A generalization to models incorporating fluc- 
tuating environmental effects such as temperature, via 
statistically mixed initial states, is still an open question 
on which work is in progress. 

The strategy is exemplified by two illustrative cases, 
the alignment and the orientation control of rigid-rotor 
molecules. We show that either complete or close to com- 
plete control is achieved after the application of a finite 
number of pulses. More precisely, the alignment / ori- 
entation control processes, aiming both at efficiency and 
maximal duration within a rotational period, require few 
kicks to get close to the target, while remaining in the 
conventionally fixed finite dimensional Hilbert subspace. 
We note that in all the dynamical mechanisms under con- 
sideration, we are interested in controlling the post-pulse 
behavior of the observables. The motivations for such 
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processes are now well established in a wide variety of ap- 
plications extending from chemical reaction dynamics to 
surface processing^ catalysis, nanoscale design, and quan- 
tum computing El [H IS • 

The paper is organized as follows. A mathematical 
setup is given in Sec. |H] for the general control strategy, 
which is then applied to the case of molecular alignment 
/ orientation. These processes are taken on a parallel 
footing in Sec. IIIII by considering either the polarizabil- 
ity or the permanent dipole moment of the molecular 
system. Section IIVI collects the results with a special 
emphasis on the robustness with respect to control pa- 
rameters and the experimental feasibility. Concluding 
remarks and prospective views are presented in Sec. E] 
Some details of the calculations and some proofs are pre- 
sented in Appendices IaI IbI IU1 andlDl 



Apart from these mathematical considerations, which 
will be used in Sec. Ill Dl the reduction of the Hilbert 
space may also lead to other physical advantages in spe- 
cific cases. Indeed, for post-pulse alignment / orienta- 
tion control, the finite dimensionality of the molecular 
rotational space, involving low momenta, leads to longer 
alignment / orientation dynamics as a result of the slower 
oscillation of (0(t)) after interaction with the field. It is 
precisely this behavior that roughly constitutes the ba- 
sis of the compromise between efficiency and long dura- 
tion of such control issues. Reciprocally, the justifica- 
tion of the dimensionality reduction is related to the fact 
that moderate perturbations, i.e., a moderate number 
of applied pulses and amplitude, can only transfer finite 
amounts of energy to the molecule, confining the system 
to a finite-dimensional subspace |l3|. 



II. METHODOLOGY 

The field-free Hamiltonian acting on the physical 
Hilbert space TL is denoted Hq. The initial pure quantum 
state of interest \tp n (t — 0)), is taken as one of the eigen- 
states of H corresponding to the eigenvalue E n . One 
considers an observable O which, through the dynamical 
behavior of its expectation value (0(t)), describes the 
physical property to be controlled by the external field. 
O is a self-adjoint operator which does not commute with 
Hq. The aim of the control is to find an evolution opera- 
tor U(t), within a class of experimentally realizable pro- 
cesses, such that the time evolution of the initial state un- 
der the action of U (t) leads to an average (ijj n (t)\0\ip n (t)) 
that is maximized or minimized. This goal is achieved in 
three steps. 

A. Finite dimensional subspace 

We first select a finite A^-dimensional subspace TiS N > 
of the physical Hilbert space. For instance, in the case of 
alignment / orientation such a subspace is generated by 
the first N eigenstates \n) of H with eigenvalues E n [l]|. 
The mathematical advantage of this reduction is twofold 
when considering the reduced operator 

q(N) = p (N) 0p (N) ^ 

pW being the projector on the subspace H^ N '. The first 
advantage, of kinematical nature [iJInj], is related to the 
fact that £>W has now an upper and lower bounded 
discrete spectrum, as opposed to the possibly continuous 
or unbounded spectrum of O. The second is that the 
controllability of the system can be completely analyzed 
(see Sec. lIIDl for details), whereas such results are limited 
for the infinite-dimensional case [22|. 

Finally, we assume that all the fundamental frequen- 
cies of Hq in 1HS N ' are commensurate, i.e., integer mul- 
tiples of a fundamental frequency lj, which implies that 
the motion is periodic in time, with period T = 2ir/u. 



B. Controllability 

The control is exerted through the application of a 
time-dependent external field described by the Hamil- 
tonian 

H(t) =H + V(t), (2) 

where V(t) is the molecule-field interaction potential. 
The general structure of V{t) = V(£, £ (t)) involves the 
field vector £ (t) (amplitude and polarization direction) 
together with molecular vectorial characteristics £ that 
may, in some cases, be field-induced (an example of such a 
behavior, through the polarizability interaction, is given 
in Sec. IHIjl . Moreover, in order to apply the results of 
control theory, we also assume that the system is control- 
linear, which means here that V can be written in the 
form 

V(t) = -v{t)Hi, (3) 

where v is a real control function (that can be chosen 
at will) and Hi the interaction operator. If a system of 
the form © is completely controllable [jllll^ then there 
exists a function v(t) such that the corresponding prop- 
agator U (t) brings the initial state to the target state. A 
simple example of a non-controllable system is a decou- 
pled one. In this case, the matrix representation of H is 
block-diagonal and the Hilbert space can be decomposed 
as a direct sum of at least two orthogonal subspaces. 
However, for non-decoupled systems, the question of dy- 
namical realizability is more difficult and an investigation 
of the structure of the Lie algebra is needed. A Lie alge- 
bra of matrices is a subspace of the vector space (over the 
field of reals) of N x N matrices (with complex entries) 
which is stable under the commutation operation |llj |. 
Two cases of particular relevance here are the Lie alge- 
bra of iV x JV skew-Hermitian matrices denoted by u(N) 
and whose dimension is A^ 2 , and the Lie algebra su(N) 
of JVxJV zero-trace skew Hermitian matrices, whose di- 
mension is A^ 2 — 1. We recall 0, ^| that a necessary 



3 



and sufficient condition for the complete controllability 
of the Hamiltonian system defined by H [Eq. |J2J] is that 
the Lie algebra generated by iHq and iHi be u(N). In 
the rest of this section, we will assume the dynamical 
realizability of the target states. 



C. Target state 

The state |x ) which maximizes (or minimizes) Q( N > 
in the subspace TiS N ' is nothing but the eigenstate of 
£)( N ' corresponding to its upper (or lower) eigenvalue. 
This provides a clear description of the target state in 
terms of a linear combination of states \n) in direct 
relation with the observable to control. If N is large 
enough, the upper (or lower) eigenvalue of 0^ will be 
close to the corresponding bound of O (if this operator 
is bounded) and an efficient control could be achieved 
(assuming complete controllability). Once again, it is 
worth noting that large N allows for higher values of 
\{OW)l but leads to higher energy transfer from the 
external field to the molecular system and consequently 
results into shorter durations as the post-pulse average 
(0(t)) is done over a superposition involving many more 
eigenstates of H . 



D. Control strategy 

We shall now present in detail two control schemes. 
The first one, which was first proposed in |23( and fur- 
ther developed in p|, consists in applying pulses each 
time (OW) reaches its maximum value within a period 
T of the field-free dynamics. An individual pulse per- 
turbs the molecule by an effective evolution operator U A , 
where A is a real parameter related to the field ampli- 
tude: Eq. (7J gives, for instance, the expression of A in 
the sudden approximation. The system remains in the 
subspace HS N \ up to negligible corrections, if the field 
amplitude A is small enough. Between two successive 
pulses, the molecule evolves following its field-free mo- 
tion, governed by [/# (Ai) = exp(— iHoAt). An impor- 
tant hypothesis of our strategy is that the system is to 
be perturbed according to a unitary operator U A , which 
commutes with 0^ N \ such that its application does not 
alter (flW) = {U^O^U A ) . Furthermore, the optimal 

target state |x^) is an eigenfunction of both £)W and 
U A . 

We now consider that pulses are applied at times 
ti when (OW(t)) reaches its global maximum Oi — 
(0^ N '(ti)). The sequence of Oi's is increasing (possi- 
bly constant) but bounded and is therefore convergent. 
Due to the hypothesis on the fundamental frequencies 
(Sec. EH) and the fact that = {U^O^U A ), 

(t)) is a periodic, continuous and differentiable (un- 
der free-evolution) function. It is then easy to verify 
that the limit of the preceding sequence is a fixed point 



Oi = Oi + i, corresponding to a wave function \ipf) (or 
a family of wave functions) such that (V>/|0( Ar )|-0/) — 
{ , if)f\U'^ 1 0^ N 'U A \'if)f) = Oi is a global maximum within 
a period T for any value of the parameter A. 

The description of this control strategy raises several 
questions which are not completely resolved. The first 
one is the set of limits (corresponding to the fixed points) 
of the process. It can readily be shown that this set con- 
tains some of the eigenvectors of the operator 0^ N \ but 
we can look for the conditions on Hq, U a and 0^ N ' such 
that these wave functions be the unique fixed points of 
the strategy. Appendix 1X1 provides a sufficient condition 
on these operators. The second question deals with the 
determination of the domain of attraction of each fixed 
point and, more precisely, the domain of the target state. 
The domain of attraction of |?/>/) can be defined as the 
set of initial states such that, for some choice of the 
parameters Ai (the pulses being labeled by i), the pro- 
cess, starting from converges to I?/ 1 /)- The last open 
question concerns the reduction to a finite dimensional 
space of the dynamics. In this paper, the fact that the 
dynamics resides within the subspace H^' is only nu- 
merically justified (see Sec. IIIII for details). Work is in 
progress on these open questions. 

The second control scheme that can be used is a strat- 
egy similar in its spirit to the first one. Up to now, we 
have searched for a sequence of perturbations that max- 
imize the average value of the observable (OW(t)), An- 
other relevant optimization would be to get closer to the 
target state \x^ N ^)- This can be done in practice by ap- 
plying pulses each time the modulus of the projection 
on the target state, |(x ' |t/; n (t))|, reaches a maximum. 
The relative merits of the two control schemes will be 
discussed in Sec. IIVI 

To conclude, we mention an alternative that can fa- 
vor the convergence of the control scheme. The general 
strategy is based on global maxima of the field-free dy- 
namics, but local maxima can also be used. This will be 
illustrated numerically for the control of molecular align- 
ment in Sec. IIVI where it will be shown that a better 
convergence is obtained with local maxima. 



E. The sudden approximation 

Depending on the duration t of the radiative interac- 
tion V(t) as compared to molecular characteristic times 
T, we may be in the adiabatic (r 3> T) or in the 
sudden (r <C T) limits. Long laser pulses have been 
shown to yield very efficient adiabatic molecular control 
pil I25L 12a. |27| , but with the caveat of it disappearing 
after the pulse, as the molecule returns to its initial state 
when the turn off is adiabatic. As we are rather con- 
cerned with the post-pulse behavior of (0(t)}, we here- 
after use short pulses, assuming the sudden limit. Fol- 
lowing the derivation of the lowest order impulsive ap- 
proximation given in Refs. [28l 12^. the time-dependent 
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Schrodinger equation, written in the interaction repre- 
sentation as 

ihjfiMt) = - [e iH ° t / n V(t)e- iH ° t / n ] Mt), (4) 

with "ij} n (t) — e lHot ^ h ip n (t), is integrated on the total 
pulse duration t to yield 



h Jo 



■Hot/H v ^ e -iH t/h\ ^ n (t) dt+ ^ n (t = 0). 

(5) 

At the lowest order in t/T, one can substitute e lHot ^ h 
with 1 in Eq. © and, using V(t) = -v(t)Hx [Eq. ©], 
the integral equation (J5J can be solved as 



i> n (r) = exp(iAHi)i> n (t = 0), 



with 



A = - 



v(t)dt. 



(6) 



(7) 



Finally, returning to the Schrodinger representation, one 
obtains 

Mt>r) = U Ho (t)U A Mt = 0), (8) 
where the evolution operators are given by 

U Ho (t) = exp(-iH t/h) (9) 

and 

U A = exp(iAHi). (10) 

We are now in a position to examine the conditions 
(choice of the parameters N, A, and tj's) under which 
the alluded strategies work in the control of alignment / 
orientation processes. 



III. CONTROL OF THE MOLECULAR 
ALIGNMENT / ORIENTATION 

Laser-induced alignment and orientation are impor- 
tant issues in controlling molecular dynamics 0- For 
a diatomic molecule driven by a linearly polarized laser 
field, alignment means an increased probability distribu- 
tion along the polarization axis whereas orientation re- 
quires in addition the same (or opposite) direction as the 
polarization vector. In a variety of applications, extend- 
ing from chemical reaction dynamics to surface process- 
ing, catalysis and nanoscale design 0, 0, 0, 0] , notice- 
able orientation that persists after the end of the pulse 
is of special importance. This is why the identification 
of target states which fulfill these two requirements (i.e., 
efficiency and persistence) is the basis of the construction 
presented in Ref. 0. 

The time evolution of the molecule, hereafter described 
in a 3D rigid rotor approximation, interacting with a 



linearly-polarized field is governed by the time-dependent 
Schrodinger equation (TDSE), expressed in atomic units, 



BJ 2 - $(£) ■ £{t) ^(M;t), (11) 



where J is the angular momentum operator, B the ro- 
tational constant, fl(£) the total dipole moment and £ 
the electric field. 9 denotes the polar angle between the 
molecular axis and the polarization direction of the ap- 
plied field. The motion related to the azimuthal angle 
4> can be separated due to cylindrical symmetry. The 
field-induced dipole is expanded in powers of £ (t), 



fio + a ■ £(t) + 



(12) 



where the two first terms (to which our approximation 
is limited) are respectively the permanent dipole fig and 
the polarizability tensor a. We notice that only very 
intense lasers may require the inclusion in Eq. 112J1 of 
higher order terms. The electric field is written as 



£(t) = f(t)cos(ut)e, 



(13) 



e being the unit polarization vector, oj the carrier wave 
frequency and f(t) the laser pulse envelope. In a high- 
frequency regime, which is often referred to in models 
dealing with alignment 2j, the permanent dipole in- 
teraction, i.e., fig • £(t), averages to zero due to the 
term cos(ui), whereas the polarizability interaction cor- 
responding to a ■ £^(t), involves a cos 2 (ujt) term which 
contributes an average 1/2. Finally, in a high-frequency 
approximation [SEfl' tne TDSE [Eq. (JTTJ] is written 



1^(9, fot) = 



BJ 2 - -f{t) (Aacos 2 



aj_ 



<p(e,4>; t), 

(14) 



where Aa = a\\ —a± is the difference between the parallel 
or|| and perpendicular a±_ components of the polarizabil- 
ity tensor. While excellent alignment can be obtained 
through adiabatic transport on so-called pendular states 
|26j | resulting from field dressing, only sudden pulses offer 
the possibility of alignment that persists after the field is 
over 0. As opposed to alignment, orientation requires 
spatial symmetry breaking, and therefore cannot be ob- 
tained from a radiative interaction in cos 2 9 such as the 
one depicted in Eq. I|14|l resulting from polarizability. It 
has recently been shown that very short pulses combining 
a frequency u> and its second harmonic 2u (in resonance 
with a vibrational transition), excite, through the com- 
bined effect of /iq and a, a mixture of even and odd ro- 
tational levels and have the ability to produce post-pulse 
orientation |3jJ. But even more efficient orientation is 
obtained using half-cycle pulses (HCPs), that, through 
their highly asymmetrical shape, induce a very sudden 
momentum transfer to the molecule which orients under 
such a kick after the field is off jH, H3 • It is worth not- 
ing that both the (lu + 2lu) and the kick mechanisms have 
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received a confirmation from optimal control schemes [7J . 
In a moderate field approximation where the polarizabil- 
ity interaction is neglected as compared to the permanent 
dipole one, the TDSE which governs the orientation pro- 
cess can be written as 

*J^(M;t)= [BJ 2 - f(t)vo cos 8] ip(6,<P;t). (15) 

As the two processes (i.e., post-pulse alignment / orien- 
tation) require short duration pulses, from now on, we 
assume a sudden approximation based on the shortness 
of the pulse duration t as compared to the molecular ro- 
tational period T rot = tt/B. For relatively low j (where 
j labels the eigenstates of J 2 ), this amounts to the def- 
inition of a dimcnsionlcss, small perturbative parameter 
e = tB. This definition, together with a rescaling of time 
s = t/r (such that s S [0,1] during the pulse) leads to 
a TDSE that can be treated by time-dependent unitary 
perturbation theory [33T.l34| . 

i^(0, 0; s) = [Ho - Va-As)} t/>(0, & a), (16) 

with Hq = eJ 2 and V" a;0 (s) the radiative interaction de- 
pending on the process to be described. More precisely, 
we have for alignment 

V a (s) =E 2 a (s) cos 2 + F 2 (s), (17) 

and for orientation 

V (s)=E (s) cos 6. (18) 

The time-dependent functions are defined as 

E 2 a {s) = Aarf 2 (T S )/2, (19) 
F 2 (s) = a ± rf 2 (r S )/2, (20) 
E (s) = W)t/(ts). (21) 

Returning now to the successive steps of our approach 
as detailed in the mathematical setup (Sec.[n]), we have 
to define an initial pure quantum state and the observ- 
ables adequately describing both processes. The initial 
state of the molecule is taken as the ground state of the 
rigid rotor, 

^„=o(t = 0)) = |i = 0,m = 0). (22) 

We remark that any other rotational state of the rotor, 
i.e., \j > m,m) could not be taken, by observing that 
the projection quantum number m is invariant upon the 
application of linearly polarized pulses. The two observ- 
ables we consider are O a — cos 2 9 for alignment and 
= cos 9 for orientation. The goal of the field driven 
molecular alignment (or orientation) is then to maximize 
(or to minimize, depending the choice of the orientation) 
for the largest time duration, the expectation values 

(O a;o (s)) = (i>(6, 0; s)\O a;o \^(0, 0; «)>, (23) 



after the pulse is over. 

The second step consists in the reduction of the dimen- 
sionality of the infinite physical Hilbert space by consid- 
ering a finite subspace H^2 generated by the first N 
eigenvectors of J 2 , i.e., \j, m = 0) (j = 0, 1, . . . , N — 1). 
The observables are projected in this subspace according 
to 

q(N) _ p{N) Q p (N) (2A) 

with the projectors explicitly given by 

JV-l 

fi= ) o = E^ m=0 )^ m=0 i- ( 25 ) 

3=0 

It is clear from Eq. I|24|) that oi^} has, as opposed to 
cos 2 9 or cos 9, a purely discrete spectrum. The target 
states Ixlfo ) which maximize the alignment (or the orien- 
tation) in the subspace Tl^l are simply the eigenstates 
of Oafo with the highest eigenvalue. Their calculation 
in terms of the appropriate expansion coefficients c°'° in 
the \ j, m = 0) basis, 

JV-l 

lxif o ) )=E c rii' TO = )' (26) 

involves the diagonalization of the [N x N) matrices of 

Oo-c? [Eq. Q24JI] written in the same basis. For instance, 
in the simplest case of orientation, a tridiagonal matrix of 
the operator Oi* results (when using the approximation 
(j, m = 0| cos9\j ± 1, m = 0) = 1/2, valid for j m) in 
the analytical expression 

i^^riH^) b ' ,ro=o) ' 

. < 27 > 

with the corresponding maximal orientation found in this 
subspace 

( x W|cos%f))^cos( ] ^ I ). (28) 

Having determined the kinematical constraints, we can 
now analyze the dynamical realizability of these target 
states. This kind of Hamiltonian (involving a dipole in- 
teraction with the control field) has already been inves- 
tigated by S. G. Schirmer and co-workers [lj and it was 
shown that the Lie algebra gene rated by iH^ and iO^J 
is u(N) (Theorem 1 of Ref. |12|). The system is therefore 
completely controllable. 

The next step consists in applying our control strategy. 
For doing so, we determine the evolution operator for 
individual pulses. This is done in the impulsive sudden 
limit |28l l29l |. The interaction term of the alignment 
process V a (s) [Eq. (|17fl ] involves the square of the electric 



6 



field amplitude. The time evolution can be expressed in 
terms of a free evolution combined with 



where 



U a = exp [iA a cos 2 9] , 



A a = [ E 2 a (s)ds 
Jo 



(29) 



(30) 



is the total pulse-square area. We note that we have not 
taken into account the ^-independent contribution F 2 (s) 
to V a (s) as it would only contribute a pure phase factor. 
For the case of orientation, we get as the instantaneous 
evolution operator describing the kick 



with 



U = exp [iA cos 9] , 



A = ( E {s)ds 
Jo 



(31) 



(32) 



the total pulse area. Between pulses, the molecule evolves 
under the effect of its field-free rotation, described by the 
evolution operator 

Uh = CX P [— «e J 2 s] . (33) 
One can easily check the necessary commutation relations 



{N) U 



0. 



(34) 



implying that the target states |xi^c?) are eigenfunctions 

of both oifo* and U a - . Note that all quantities (s, e, 
E 2 , E D ) are dimensionless, and hence do not depend on 
a specific molecular system. Times can be expressed in 
terms of fractions of the molecular rotational period T rot . 
Moreover, the molecule-laser interaction characteristics 
(E 2 or E ) combine the field strength £ with molecular 
constants a or fiQ, in such a way that the same values 
may be considered for different molecules at other, ade- 
quately chosen laser intensities. The final step for the 
control scheme consists in applying the instantaneous 
perturbations at specific times. More precisely, as has 
been discussed in the mathematical setup (Seem, two 
distinct strategies are adopted for the determination of 
these times. The first one (SI) consists in taking the se- 
ries of Si when (0^ r J)(s) reaches its successive maxima 
(global or local) during the free evolution following a kick. 
The second one (S2) consists in taking the s,'s as times 
when the modulus of the projection of the wave function 

on the target states \(ip\(s)\xa?o)\ reaches a maximum 
during the free evolution. The numerical results obtained 
by both strategies are presented in Sec. IIVI 

We now discuss the relative merits of the two schemes 
from a theoretical point of view. For the strategy SI, the 



slope after the pulse can be explicitly determined in the 
infinite-dimensional shown in Appendix O 

d , 



ds 



(COS ( 



= 2eA(l - (cos 2 9}), 



ii+O 



for orientation, and 



ds 



2Ae(sin 2 29), 



(35) 



(36) 



Si+O 



for alignment. As (cos 2 9) < 1 and (sin 2 29) > 0, one de- 
duces that, in both cases, the slope takes a non-zero value 
after the application of the kick, whatever the pulse area 
A of the field. This point proves that the control scheme 
can be iterated for the alignment / orientation processes. 
In the finite dimensional case, the corresponding expres- 
sion for orientation derived in Appendix O is slightly dif- 
ferent. We obtain for cos^ 9 

Wfl) = 2eA(l-((cos^9) 2 )+B)- 

Si+O 



ds 



(cos 



0(A 2 ), 
(37) 



B being a boundary term which is of the form 

B = 2N +1 M ' (38) 

where the a n 's are the coefficients of the wave function 
in the basis |0 < j < N — l,m = 0). Here, we notice 
that we cannot conclude as above and arguments of di- 
mensionality have to be used. In this way, Appendices 
1X1 and iBl show that the unique fixed point of strategy S2 
is the target state |x^) whereas scheme SI can possess 
several fixed points as the eigenfunctions of the opera- 
tor 0^ N \ However, the strategy SI has the advantage 
that the target state can be selected during the dynami- 
cal process by appropriately choosing the intensities and 
the number of laser pulses. This is not the case for the 
scheme S2 because the target state has to be selected a 
priori. 



IV. RESULTS 

The results are presented in a general way, transfer- 
able to any particular molecule (characterized by its ro- 
tational period, its permanent dipole and its polarizabil- 
ity) interacting with any pulsed field (characterized by 
its moderate amplitude and its short duration as com- 
pared to the molecular rotational period). The relevant 
parameters which gather both molecule and field con- 
stants are taken to be A a — 1.5, A — 1, and e = 0.03. 
They are compatible with any molecule, provided that 
the field amplitude is adequately chosen. For LiCl, for 
instance, this amounts to a pulse duration of about 0.3 
ps and a field amplitude of 1.5 x 10 5 Vcm -1 |2£j in the 
case of the orientation. Time is indicated in fractions of 
the molecular rotational period. The presentation of the 
results follows the different steps of the control strategy. 
Analytical estimations of the dynamical parameters are 
reported in Appendix iDl 
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be done at that step. In order to keep an alignment / 
orientation duration of the order of one tenth of the ro- 
tational period (amounting to durations exceeding 10 ps 
for heavy diatomic molecules such as Nal), N has to be 
limited to 5, as can be observed from Figs.^ an d0 This 
may seem rather limiting in view of the moderate molec- 
ular rotational excitation, but even such a low dimension 
(JV = 5) turns out to be sufficient for very efficient align- 
ment (about 0.85) or orientation (about 0.9). This fixes 
completely the target states (approximately given by Eq. 
(|T7Jl in the case of orientation) for both processes. 



B. Dynamics controlled by the two strategies 



FIG. 1: Maximum alignment efficiency (crosses) and associ- 
ated duration (open circles) as a function of N, the dimension 
of the rotationally excited subspace 7i^2 (see text). The 
solid and dashed lines are just to guide the eye. 




23456789 10 11 

N 



FIG. 2: Same as FigQ but for orientation. 



A. The finite dimensional subspace 

Relevant information for the characterization of align- 
ment and orientation in the finite subspace Tx^Iq as a 
function of its dimensionality N are displayed in Figs.^ 
and [21 The maximum efficiency (approximately given 
by Eq. I|28(l for the case of orientation) that can ideally be 

expected for a process that stays confined within TC^Iq 

is given as (xo |Oa-c? |Xo )■ The niaximum duration, 
measuring the relative duration of the alignment / ori- 
entation processes over which (oi%\s)} remains larger 
than 0.5 during the field-free evolution, is indicated in 
terms of a fraction of the rotational period as At/T rot . 
The information displayed is crucial for the choice of the 
dimensionality N of the reduced subspace: a larger N 
is more suitable for a better efficiency of alignment / 
orientation, but leads to a shorter duration. The com- 
promise between maximum efficiency and duration is to 



The two strategies, noted 51 and 52, can be exploited 
in the same way, the time delays between the successive 
pulses being defined following one of the two processes 
(maxima of (Oa-o) for 51, or maxima of | (xa$>\ip\(t)) \ 2 
for 52). We look for the dynamical behavior of the 
two quantities, namely the expectation value of the ob- 
servables (O^l) and the projections on the target states 

\(Xcr)>\'>P\(t))\ 2 ■ All results are displayed with the first 
pulse taken as the origin of time and extend at least one 
rotational period T rot after the last pulse. This helps to 
show the complete field-free behavior of the dynamics, 
which repeats periodically leading to revival structures. 
Figures |3 and 01 illustrate the results of applying strategy 
51 to the alignment and orientation processes, respec- 
tively. Noticeably different behavior is obtained for the 
two processes: clearly, a smaller number of pulses is nec- 
essary to control alignment (i.e., 6 pulses for alignment 
instead of 15 for orientation). The first pulse already 
provides an alignment efficiency of 0.75, whereas it only 
yields 0.5 orientation efficiency. Conversely, the over- 
all dynamics (conditioning the delays between successive 
pulses) is much more oscillatory with complex structures 
in the case of alignment as compared to orientation where 
it progressively increases up to 0.89, which is almost the 
optimal limit as found from Fig. [21 In ah cases, a com- 
parison of the observables and the projections calculated 
using the exact wave function ip\{s) with the one prop- 
agated in the subspace 7i. m=0 supports the claim that 

(5) 

the rotational dynamics actually resides within H. a . In 
both cases, panels (b) shows the way the wave function 

tpx(s) gets close to the target state xi^o, illustrating the 
successful outcome of the control strategy with a coher- 
ent choice of parameters A, N, and the number of pulses 
for appropriately describing the dynamics. Finally, the 
post-pulse dynamics leads to results particularly remark- 
able with respect to previous proposals. An efficiency of 
about 0.85 with a duration of (l/10)th of the rotational 
period is achieved for alignment, with even better results 
for orientation, i.e., 0.89 efficiency and (2/10)th of the 
rotational period duration (which corresponds to about 
2 ps for a light molecule like LiCl and 20 ps for a heavy 
one like Nal). 
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Q l , , 1 1 , L 

oi < b ) 



0.3 0.6 0.9 1.2 1.5 1.8 

t/T , 
rot 



t/T , 
rot 

FIG. 3: Alignment dynamics during and after the train of 
pulses determined by strategy SI of perturbing at the global 
maxima of (cos 2 9) : panel (a) for {ip\{s)\ cos 2 9\ip\(s)) and 
panel (b) for |(Xa 5 '|' ( /'A(s))| 2 - The solid line corresponds to the 
exactly propagated wave function ip\(s) and the dashed line 
to the propagation of ip\(s) in the subspace "H^L . The train 
of pulses is displayed on panel (b) and the optimal alignment 
in Ti^_ Q i s indicated by the horizontal line on panel (a). 



We next notice a series of local maxima in Fig. El (for 
instance, just after the first kick), which could be used 
for the control strategy. This point is illustrated in Fig. El 
for strategy 51. It turns out that excellent results are 
obtained for the alignment process using in the control 
scheme the local maxima (the next one after the appli- 
cation of the pulse) . More unexpectedly, this leads to a 
better result than the use of the global maxima, as opti- 
mal alignment is almost achieved after only 4 kicks (while 
6 are used in Fig.EJ). In this example, it is clear that the 
best choice for the convergence of the process is to kick 
at a local maximum. However, we did not manage to 
establish general rules for this degree of freedom and it 
thus seems that numerical tests have to be undertaken 
for each practical case. 

The outcome of the calculation when applying strategy 



FIG. 4: Same as Fig. |3] but for the orientation dynamics. 

52 to the orientation, taken as an example, is displayed in 
Fig- El The results obtained are very similar to those us- 
ing strategy 51, with a smaller number of pulses (viz., 9) 
but leading to slightly less efficient orientation. This sug- 
gests a great similarity between the two strategies from a 
numerical point of view, as can be seen when comparing 
Figs. Hand El 



C. Robustness 

The robustness of the overall strategy has in princi- 
ple to be checked against two variables: the time delays 
between successive pulses and the total energy delivered 
by individual pulses. Figure illustrates the robustness 
against the most sensitive parameter, i.e., the time de- 
lays Si, where (Oi N \si)) is maximum. This is given for 
the orientation process only, by considering again the two 
strategies 51 and 52. Strategy 52 seems here more ro- 
bust, but this is only due to the smaller number of pulses 
which are necessary for reaching the target state (9 for 
52, instead of 15 for 51). 

As for experimental feasibility, the accuracy of the light 
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0.3 0.6 0.9 1.2 



t/T , 
rot 

FIG. 5: Same as Fig. [3] but local maxima (the first maximum 
after the pulse) have been used. 

paths followed by successive pulses is of the order of 
0.2 /im, corresponding to 0.6 fs. This is, for instance for a 
light diatomic molecule like LiCl, less than l/10th of the 
inaccuracies considered when varying the time delays, in 
Fig. {7\ In other words, both strategies are very robust 
with respect to such inaccuracies affecting the molecule- 
pulse interaction terms. On the contrary, the measure 
of the total energy per pulse is affected by a more se- 
vere inaccuracy. This is checked in Fig. |SJ by varying 
A by 10%, for both strategies 51 and 5*2. Remarkable 
robustness is achieved, advocating for an experimental 
feasibility of the control scheme: a few short pulses (say 
5 to 6), with a moderate total energy determined within 
10% of accuracy, applied with time delays of the order of 
1 ps with an accuracy of 1 fs, leads to an excellent control 
of the orientation for any molecular system. 



V. CONCLUSIONS 

In this work, we have developed a laser control strat- 
egy for a molecular system initially in a pure quantum 



1 




state, aiming at the maximalization (or minimalization) 
of an observable O given by an operator, not commuting 
with the molecular field-free Hamiltonian. The strategy 
consists of four steps: reduction of the dimensionality of 
the original physical Hilbert space; definition of a tar- 
get state as an eigenvector of the reduced observable in 
the finite subspace; analysis of the controllability of the 
system; application of a train of sudden perturbations U, 
commuting with O, at times where the expectation value 
of O reaches a maximum during its field- free evolution. A 
process similar in its spirit consists in applying the pulses 
at times when the projection of the time-evolved wave 
function on the target states reaches maximum values. 
This latter strategy yields comparable numerical results 
even if, from a theoretical point of view, its set of fixed 
points (which corresponds only to the target state) can 
be smaller than (or equal to) the one of the first process. 
Moreover, it is worth noting that such general strategies 
are not only independent of the particular molecule, but 
are also transferable to different control issues and differ- 
ent observables. Finally, we have checked the robustness 
with respect to inaccuracies in the time delays between 
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FIG. 7: Robustness of the orientation dynamics when shift- 
ing the time delays between pulses systematically by ±0.1% 
of the rotational period (the dashed line corresponds to over- 
estimates by +0.1% whereas the dotted-dashed one to un- 
derestimates by —0.1%, the solid line corresponding to the 
calculation of Fig. 2J. The strategy SI is displayed on panel 
(a) and S2 on panel (b). 



FIG. 8: Robustness of the orientation dynamics varying the 
pulse energy by 10%. The full line corresponds to A = 1, i.e., 
the calculation of Fig. [I] the dotted and dotted-dashed lines 
are for A = 1.1 and A — 0.9, respectively. The strategy 51 is 
displayed on panel (a) and 52 on panel (b). 



pulses or the total pulse energies, and found it to be re- 
markably good, as long as moderate intensity and short 
radiative interactions are considered. 

Control strategies based on pure states are however 
subject to important modifications when dealing with 
temperature effects. Orientation, for instance, drastically 
decreases for increasing temperatures j33, HE HE l37j . 
This is basically due to the fact that, at non-zero temper- 
ature, the initial state is in general not a pure state, but 
a mixed one, i.e., a superposition of a statistical ensem- 
ble of rotational states with m ^ 0, which tend to mis- 
align the molecule. The dynamics of such mixed states 
is described by a density operator p(t) (instead of a wave 
function) evolving according to the von Neuman equation 
(instead of the Schrodinger equation) . Following the di- 
mensionality reduction of the Hilbert space, the control 
strategy looks now for the optimal density operator 



which maximizes the observable Q( N > in the same sub- 
space. This amounts to maximizing Tr{p( N >0( N >], with 
the commutation relation [p^ N ',0^ N '] = 0. It has been 
shown [E E3 that such a goal can be achieved by an uni- 
tary operator U that appropriately orders the common 
eigenvectors of and 0^ N \ such that the one corre- 
sponding to the highest eigenvalue of p( N ^ be associated 
with the one corresponding to the highest eigenvalue of 
0^ N \ and so forth. An even more sophisticated gener- 
alization of the control strategy is to be looked for when 
the particular system under consideration evolves while 
interacting with a physical environment describing colli- 
sional processes in a gas, or friction in a liquid |38l l39j. 
The strategies that have to be worked out would involve 
non-unitary perturbations U and a special attention has 
to be paid to decoherence [iol l4l| . Work in these direc- 
tions is in progress. 
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APPENDIX A: ANALYSIS OF THE SET OF 
FIXED POINTS: THE GENERAL CASE 

The goal of this appendix is to determine the set of 
fixed points of the control strategy SI. Indeed, as men- 
tioned in Sec. [n] a crucial question for the convergence 
of the control scheme is the set T of its fixed points. In 
general, we are not able to exactly determine T but only 
a larger set S. The latter contains the wave functions 
\tp s ) which fulfill the requirements 



and 



(M#o,0 (JV) ]|^) =° ( A1 ) 



(i> s \Uf[H ,OW]U A \Tp s )=0 ) (A2) 



for all values of the real parameter A. As mentioned 
above, we can show that T is a subset of S. For the 
purpose of proving this result, we observe that just before 
the pulse at ti — one has 



dt 



0. 



whereas just after the pulse at ti + 



ti+O 



= 1 {u- a 1 [h M n) ]u a ) 



= i ([U A 1 H U A ,O^}) 



ti+O 



ti+O 



(A3) 



(A4) 
(A5) 



We can conclude by arguing that if there exists a value of 
the parameter A such that the slope undergoes a change 
from zero [Eq. (|A3I) ] to a non-zero finite value [Eq. (|A4I) ]. 
then (C (Ar) (t)) will reach within the period T a maximum 
strictly larger than the one obtained prior the application 
of the pulse. 

We are now in a position to state the following theo- 
rem: S = \\x^)} (the |x^)' s being the eigenstates of 
£>W) if and only if the dimension of the vector space 
V generated by [Hq,0^] and the elements of the form 
Uj^HaiOWpfr where A G M, is equal to N(N - 1). 
The proof goes as follows. 

We first recall that the complete controllability of the 
system is assumed. This condition is equivalent to re- 
quiring that the Lie algebra C generated by iHq and iHi 
is u{N) or that the dimension of this vector space is N 2 . 
Moreover, {/^[Hq, 0^]U a being a skew-hermitian op- 
erator, it can be shown that V is a subspace of C. The 
starting point for the following computations will be the 



where e mn = \x {m} )(x (n) \, I < n < N — 1 and n < 
m < N . We denote D the dimension of V. This vector 
space is generated by the set Vk (k — 1, . . . , D), which 
can be written as a linear combination of elements of 
the basis given by Eq. l|A6|l. Expanding [H ,O^] and 
f/^" 1 [i?o, 0^]U A with respect to the previous set and 
taking the mean value of these operators over a wave 
function \ip s ) S 5, we obtain a system of linear equations. 
Inverting the system, one gets D equations of the form 



(n)\ 



G S and ( X (n) \d n , 



\X 



0. 

(«)\ 



(A7) 



1 for any 



As ... 

1, . . . , N, it is clear that V does not contain any linear 
combination of the d„ n 's, which shows that the dimen- 
sion of V is at most N(N — 1). Hence, if we assume 
that, for instance, x mn G V and y mn G V, we obtain by 



introducing the wave function 
Eq. EH 



^m. ~l~ Cr> 



T N c 



IX 



(n)\ 



into 



= 0, 
= 0. 



(A8) 



From Eqs. (|A8|) . one deduces that c m = or c n = 0. By 
repeating this procedure for each couple (m, ri), it can be 
shown that the states |x^), for which for some no 



(A9) 



are the unique elements of S if and only if the dimen- 
sion of V is N(N — 1). Numerical calculations of this 
dimension and complementary results are presented in 
Appendix [B] Finally, a conclusion regarding this theo- 
rem is that if the dimension of V is N(N — 1), then the 
set T of fixed points of the control strategy is a subset 
of {|x^}}> the set of eigenvectors of the impulsive prop- 
agator U A , which are also the ones of as the two 
operators commute. 



APPENDIX B: ANALYSIS OF THE SET OF 
FIXED POINTS IN THE SUDDEN 
APPROXIMATION 

This appendix particularizes the general arguments of 
the Appendix to the sudden approximation. We first 
consider the strategy 51. In this case, we have to de- 
termine the dimension of the vector space V defined in 
Appendix^] In the sudden approximation, we recall that 
the evolution operator U A can be written as [Eq. (fTT3f> ] 



AAHi 



(Bl) 
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and V is then generated by the elements of the form 



-iAHj 



[H ,OW], 



iAHr 



where A £ R. For the sake of 
clarity, we now consider that Hi = (this assump- 

tion is fulfilled for the alignment / orientation processes). 
Using the Campbell-Hausdorf formula simple al- 

gebra shows that V is also spanned by the operators 
ad n (i? ,O (A,) ) (n > 1). The notation ad n (B,C) for the 
operators B and C is defined by the recurrence formula 



ad n {B,C) 



B n = 

[a^-H^C^C] n> 1 



(B2) 



From a practical point of view, owing to the complexity 
of analytical calculations, a numerical algorithm, which 
closely follows the one for complete controllability |43j . 
can be used to determine the dimension of V. The first 
step consists in rewriting the matrix of each operator 
ad n (Ho,0'' N ') (n > 1) as an JV 2 column vector obtained 
by concatenating its columns. We next construct a new 
matrix R with these vectors, each vector corresponding 
to a column of R. Finally, the dimension of V is given by 
the rank of R. 

We apply this method to the orientation process in a 
finite dimensional space. The results of these computa- 
tions are presented in Tab.[H where the notation dim(V) 
corresponds to the dimension of the vector space V. For 



TABLE I: Dimensions of the vector spaces L and V as a 
function of N for the orientation dynamics. The third column 
indicates the maximum of the dimension of V. 



N 


dim(£) 


Max. ofdim(V) 


dim(V) 


3 


9 


6 


4 


4 


16 


12 


8 


5 


25 


20 
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To conclude, we will use the following notion. A group 
of energy levels of Q^ N ' is said to be equally spaced if 
there exists a set of integers {p, q, r, s} e {1,2,..., JV} 4 
such that 



Xp - Xq = Xr - Xs- 



(B6) 



We next suppose that the spectrum of 0^ N ' is not equally 
spaced (according to the previous definition) and that the 
coefficients C mn; m^o 7^ 0. Using the fact that the left- 
hand-side of Eq. (|B5|) is a regular function of the param- 
eter A that is constant if and only if all its coefficients 
c* m c n C mn are zero, we see that the unique solutions of 
Eq. i|B5() are the eigenvectors |x ) of for which 



Cn 1 ; 



(B7) 



It is interesting to note that the two previous hypothesis 
also give practical conditions on 0^ N > and Hq in order to 
determine the set S. Indeed, for the orientation process, 
we notice that the spectrum of cos*^ 9 being equally 
spaced (the spectrum is symmetric) , the dimension of V 
is not equal to the maximum dimension JV(JV — 1). 

Let us now determine the set T of fixed points for the 
strategy S2 in the sudden approximation. In this case, 
we first clarify the definition of S. It can be shown that 
iS contains the wave functions \ip) which fulfill 



^5 



(X 



(JV)| e iAO^ 



-iAO ( 



Ho\x 



(JV)\ 



(B8) 



for any value of A, with 3 corresponding to the imaginary 
part. Eq. I|B8I) can be rewritten as 



the orientation, as dim(V) < N(N — 1), one deduces that 
S^{\x (n) )}- This point can also be understood through 
the following calculation. 

We introduce the operator C = [Hq, 0^ n '\ and we note 
C mn its matrix elements in the basis {Ix^)}) 



(x in) \c\x {m) ) 



(B3) 



For operators Hq = J 2 and = cos^ 8, it can fur- 

thermore be checked that C mn = if m = n and C mn ^ 
otherwise. Assuming that the slope is zero before and af- 
ter the pulse, we can easily show that 



E 



0. 



-"m ^ n 



o, 



(B4) 



where the c„'s are the weighting coefficients of the wave 
function just before the pulse in the basis of the eigen- 
states of Q( N \ and the Xn's the eigenvalues of 0^ N \ 
From Eqs. I|B4(I . one then deduces that 



A( Xm -Xn), 



(B5) 



5>[ 



-iA{xm-XN) 



0, 



(B9) 



the hmn's being the matrix coefficients of Hq in the ba- 
sis {Ix 1 '™'')}- If the spectrum of is completely non 
degenerate, one deduces from Eq. I|B9I) that 



-iA(Xm-XN) 



= 



(B10) 



for any value of m (m ^ JV). Simple algebra shows that 

CNC* m h mN = 0. (Bll) 

If we now assume that h m N 0, one obtains that cjy = 
or c m = 0. The case cat = corresponds to the mini- 
mum of the projection \{ip\x )| 2 j which is not relevant 
here. Finally, in the other case, we can conclude that the 
unique element of S is Ix^')- We also remark that these 
two conditions are fulfilled for the orientation process. In 
this example, we have therefore shown that the limit of 
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APPENDIX C: THE MOLECULAR ALIGNMENT 
/ ORIENTATION CONTROL STRATEGY : THE 
INFINITE AND THE FINITE DIMENSIONAL 
CASES 

In this appendix, we shall focus on the convergence 
of the sequence in an infinite and a finite dimensional 
Hilbert space, as only the finite case has been treated 
in Appendices and [BJ We first consider the infinite 
dimensional case. We will explicitly derive Eqs. (|35|l and 
(|36|l . For doing so, we recall the following fundamental 
commutation relations , which will be used below: 



[J 2 



cost 



[erg , cos 9] 
[J 2 , cos 2 9] 
[a, cos 2 9] 
[cos 9(j q , cos 2 9] 



2(<j0+cos9), 
cos 2 9- 1, 
4cos6lcr + 6cos 2 6>-2, 
2 cos 3 9 - 2 cos 9, 
2 cos 4 6»-2 cos 2 9, 



(CI) 



where erg = sin 9-^. We notice that such identities are 
not fulfilled in a finite dimensional subspace. We consider 
the strategy 51 for the alignment and the orientation 
processes. We assume that the average (cos m 9) reaches 
a maximum at a time Si, 



ds 



cos 



i([eJ 2 



(C2) 



where m = 1 for the orientation and m = 2 for the 
alignment. A sudden pulse being applied at the time s,, 
we have to calculate the slope 



ds 



9 e[J 2 ,cos m 9}e 



Si+O 



(C3) 

Using Eqs. IjClfl and l|C3|) . one readily obtains for the 
orientation 



ds 



i(e[J 2 , cos6»]) + 2eA(l- (cos 2 9)). (C4) 



Si+O 



For the alignment process, we have 

^-(cos 2 9) 
ds Si +o 

which can be rewritten as 

d , n 



i (e [ J 2 , cos 2 9} ) - 4 Ae ( [cos 9a e , cos 2 9} ) , 



ds 



i(e[J 2 , cos 2 9]) + 2Ae(sm 2 29). (C6) 



From Eq. (|C2(I and from the fact that (cos 2 9) < 1 and 
(sin 2 29) < 1, one deduces that, in both cases, the slope 
undergoes a change from zero to a finite non-zero value 
when a pulse is applied, for any amplitude A and any 
maximum of the average of cos m 9. Note that the sign 
of the slope and the position of the next local maximum 
depend on the sign of A. 



We address now the case of finite dimensionality. In 
such a space, the commutation relations of Eq. (|C1|) can- 
not be used and we have to introduce the following stan- 
dard basis of su(N) (related here to the eigenstates 
of H ), 



fori — ^rm ^n+l,n+l)i 



(C7) 



where e mn = \m)(n\, 1 < n < N — 1 and n < rn < N. 
\n) is the eigenvector of H with eigenvalue E n [Eq. J2)]. 
H and the interaction term Hj can be rewritten with 
respect to this basis as 



#0 - En=l E n e 



E n+ i = n(n + 1) 



where 



and 



d n +l — 



n+1 



v/(2n+l)(2n + 3) 
Assuming that ([Hq, Hi]) = 0, we now calculate 
d 



dt 



Hj 



-iAH 



'[Ho,Hi]e 



iAH! 



(C8) 

(C9) 
(CIO) 

(Cll) 



i.e., the derivative of the function (Hi)(t) after a sudden 
pulse. 

A straightforward calculation leads to the identities 



N-l 



[Ho, Hi] — — ^nd n X n ,n+l, 



(C12) 



where /i n = E n+ \ — E n , and 



1 



-[[H , Hi], Hi] = Hj 



N-l 



( 2n + l)rf, l e„+i,T l +i 

71=1 

N 

5](2n + 3)d 2 e nn . (C13) 



(C5) From Eq. !|C10|) . one can rewrite Eq. I|C13|I as 



l -[[H ,Hi],Hi] = g? - 1 + jV2 2 ^ 1 + 1 e NN . (C14) 

After expansion according to Campbell-Hausdorf's for- 
mula, one gets for the derivative of Eq. I|C11|) : 

— (Hi) = 2A(1 - (V 2 ) \a n r) 

. {N 2 + 2N + 1)N 2 , 
, „ - t== =A (a N a N ^ x - a N _ x a N ) 



(2Af + l) 3 / 2 V2iV - 1' 
-0(A 3 ), 



(C15) 
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where the a n 's are the coefficients of the wave function 
in the basis of |n)'s. We notice that Eq. 1(J15|) is very 
similar to Eq. I|C4|) except for a few terms. If A is small 
enough, the evolution of the average (Hj) will depend on 
the boundary population \a n \ 2 . 



APPENDIX D: ANALYTICAL ESTIMATIONS OF 
DYNAMICAL PARAMETERS 

In Sec. lIVI the dynamics were investigated through nu- 
merical tests. In this appendix, we give a few analytical 
estimations of dynamical parameters such as the time 
between pulses. Nevertheless, it is noted that, owing to 
the complexity of complete analytical calculations, only 
rough estimations can be obtained. 

We consider the orientation dynamics and the strategy 
51. The time At between two pulses is assumed to be 
small enough, which means that only the last kicks of the 
sequence will be correctly described (This point is clearly 
illustrated in Fig. QJ. We now analyze the evolution of 
(cos0)(s) between the times Si (taken as for simplicity) 
and Sf , a sudden pulse being applied at Si . We have 



(cos 9)(s) 



U-iAcos0 e ieJ s CQs( 



-is J s^iAcosi 



>, (Dl) 



where s = t/r is the rescaled time. Expanding the two 
time exponentials and neglecting terms of order greater 
than 2 in s, (cos (9) can be rewritten in the form 



(cos0)(s) = (cos0)(O) 



— iA cos 9 r 7 2 ___ /ji AA cos 6\ 

( cos(7je ) 
J 2 , [J 2 ,cos0]]e" 



-ies{e "™" [J , cos 0|e 
1 
2 



__ e 2 g 2 /g-iAcosSf 72 \ j2 „„„ Q^„iA cos B\ 



(D2) 



As ([J 2 , cos0])(O) = 0, the proof of Appendix iBl shows 
that 

i£ S {e- iAcose [J 2 ,cQs6]e iAcose ) = 2eAs(l - (cos 2 0)). 

(D3) 

We next simplify the third term of Eq. (|D2(1 . Using the 
commutation relations [44| 



[J 2 ,cr e ] = 2 cos J 2 , 

[J 2 , COS 0] = (70 + COS 0, 



where oe = sin0^, one arrives at 



(D4) 



(e- lAcos9 [J 2 ,[J 2 ,cos9}}e lAcose ) = 4(cos0J 2 ) 

-I- 8«^(cr 9 cos0) 

+ 4A 3 (cos0- cos 3 0) 

+ UA(1+ (cos 2 0)). 

(D5) 



Introducing the operators 

f 5R[(a e cos0)] =-i(l + (cos 2 0)), 

1 $s[(cre cos 0)1 = i-(ag cos0 + cosOag + 2 cos 2 i 



(D6) 



where 5ff and 9 correspond to the real and the imagi- 
nary part of a complex number, respectively, Eq. I|D2|) 
becomes 

(cos0)(s) = (cos0)(O) + 2eAs{\- (cos 2 0)) 

- ie 2 s 2 [4(cos0J 2 ) - 8A5((<j e cos9)) 



4A 2 (cos - cos J 



(D7) 



Equation (|D7(1 allows us to explicitly determine the time 



As = Sf — Si 



2A(1 



e[4(cos0J 2 ) -8A3((cr e cos0)) + 4A 2 (cos - cos 3 



(D8) 



Up to this point, the calculation involves no approxima- 
tion on the area of the field A, as only the time delays 
between the pulses are assumed to be small. In order to 
highlight different dynamical behaviors, we now consider 
the two limits ^4^1 and A<^1. It can be shown that 



As 



if A > 1 and 



1 



(cos 



2eA((cos0- cos 3 0)) 



A(l 



(COS' 



0)) 



2e(cos0J 2 ) 



= k < A 



(D10) 



if A 1. fc> and k K are two constants that do not de- 
pend on A. They can be estimated by assuming that the 
wave function is close enough to the optimal state |x )■ 
The first limit corresponds to the strategy suggested in 
Ref. |23, where the area A increases with the number 
of kicks. From Eq. 1D9I) . we note that the expression of 
the focusing time |23l l45l |46| is derived via a quantum 
calculation, this time being obtained in Ref. |46( through 
classical or semi-classical analyses. 

Due to the moderate intensity of the pulses used in 
our strategy (for instance, A = 1 in Fig. 0J, it can be 
approximately described by Eq. HD10|) . We recall that 
the choice of the strength of the kicks is crucial for re- 



maining in the finite dimensional subspace Ti. 



m=0- 



This 



point is clearly illustrated in Fig. which gives a view of 
the orientation dynamics under the effect of a train of 30 
HCPs of area A = 1. For the last ten pulses, the arith- 
metic average of the time delays between successive kicks 
is of the order of At/T rot ~ 5.6 x 1CP 3 . This result can 
also be derived from the analytical formula of Eq. (ID 1 Oil . 



Using the eigenstate \x^), we obtain k < ~ 5.10 3 and 



— = -As = fc<A = 5.10- 3 , 

Trot 7T 



(Dll) 



a,s A — 1. Finally, it is noted that we can also esti- 
mate the difference A(cos0) = (cos 0)(s/) — (cos0)(si) 
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